Possible Constraints on the Duration of Inflationary 
Expansion from Quantum Stress Tensor Fluctuations 



/J] 



rcJil 



Chun-Hsien Wii!l and Kin- Wang 
Institute of Physics 
Academia Sinica 
Nankang, Taipei 11529 

S' 

O ' Taiwan 
O 
(N 

5h ■ L.H. Fore 

' Institute of Cosmology 

l> 
(N 

\ Tufts University, Medford, MA 02155 

> 

(N ■ (Dated: February 7, 2008) 

o ■ 

O 

OO 

o 

\^ \ We discuss the effect of quantum stress tensor fluctuations in deSitter spacetime upon the 



Department of Physics and Astronomy 



Abstract 



O 

> 
X 

. which depend upon the duration of the inflationary period. In the context of our model, 

the effect may be used to place upper bounds on this duration. 

PACS numbers: 98.80.Cq, 04.62. +v, 05.40.-a 



expansion of a congruence of timelike geodesies. We treat a model in which the expansion 
fluctuations begin on a given hypersurface in deSitter spacetime, and find that this effect 
tends to grow, in contrast to the situation in flat spacetime. This growth potentially leads 
to observable consequences in inflationary cosmology in the form of density perturbations 
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I. INTRODUCTION 



Quantum fluctuations of tlie stress tensor operator liave been studied in numerous 
recent papers [H, i, i, i, H, 0, 0, i, i, M, M, H Q 0, M, M, M,M,^- Fluctuations 



of the stress tensor drive passive fluctuations of the gravitational field, which are to be 
distinguished from the active fluctuations due to the quantization of the gravitational 
degrees of freedom. Stress tensor fluctuations play a crucial role in stochastic gravity, 
and their role in the early universe has been discussed by several authors. 

In the present paper, we will be concerned with the effects of quantum stress tensor 
fluctuations in inflationary cosmology. Since the pioneering paper by Guth [20*], the 
inflationary paradigm has been extensively developed and now seems to be in good 
agreement with observations j2l|, 2^- One of the successes of inflation is a natural 
solution of the horizon problem; inflationary expansion allows the entire observable 
universe today to have arisen from a region within which all parts were once in causal 
contact. The rapid expansion smooths any initial classical perturbations, and leads 
to a subsequent universe which is relatively independent of the duration of inflation, 
so long as there is inflation by at least a factor of 10^^. However, it can be shown that 
inflation could not have had an infinite duration in the past [23| . A key prediction of 
inflationary cosmology is a nearly flat spectrum of initial density perturbations, which 
arise from the intrinsic quantum fluctuations of an inflaton field. However, quantum 
stress tensor fluctuations of all quantum fields, not just the inflaton, should also 
contribute to density perturbations by means of passive metric fluctuations. Unlike 
the fluctuations of a nearly free inflaton field, quantum stress tensor fluctuations are 
expected to have a non-Gaussian probability distribution. Evidence for this comes 
from calculations in simple cases which reveal that in general the third moment is 
nonzero, and hence the probability distribution cannot be symmetric 25 



Here we will consider the effects of stress tensor fluctuations of the conformally in- 
variant scalar field and the electromagnetic field in deSitter spacetime upon geodesies 
of test particles. We will employ the Raychaudhuri equation to calculate fluctuations 
in the expansion ^ of a congruence of comoving timelike geodesies, and then use the 
results to draw inferences about density perturbations in the post-inflationary period. 
The outline of this paper is as follows: Some basic formalism will be developed in 
Sec. ini This formalism will be applied to inflationary cosmology in Sec. Illli We show 
that 6'-fluct nations build up during inflation and influence the redshifting of radiation 
after inflation. In this model, we will make no specific references to the mechanism 
by which inflation ends, but will show that the resulting density perturbations grow 
as the length of the inflationary epoch increases. This will lead to an upper bound 
on the duration of inflation. In Sec. IIVI we will examine the effects of ^-fluctuations 
on the dynamics of an inflaton field, and show that there is a further mechanism by 
which stress tensor fluctuations create density perturbations sensitive to the length 
of inflation. This model leads to a stronger bound on the duration of inflation. The 
results will be discussed in Sec. |Vl 
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II. BASIC FORMALISM 



A. The Raychaudhuri equation and the conservation law 

The key tool which we will employ for studying the effects of passive metric fluctu- 
ations in deSitter spacetime will be the Raychaudhuri equation, which can be written 
for a congruence of timehke geodesies with four- velocity as 

fjf) 1 

Here 6 = is the expansion of the congruence and A is an affine parameter. In 
addition, i?^,^ is the Ricci tensor, a^" is the shear, and u;^'^ is the vorticity of the 
congruence. The vorticity may be set to be zero, and we will assume that the shear 
is negligible. In this case, the equation reduces to 

Consider the case of a Robertson- Walker spacetime, where the metric can be written 
in terms of the scale factor a{t) as 

ds'^ ^-df + (t) {dx^ + dy^ + dz^) . (3) 

In the case of comoving geodesies which remain at rest in these coordinates, — 5^ 
and the expansion is given by 

9 = 00^3-, (4) 
a 

where d — da/dt. With a perfect fluid source, the matter stress tensor is 

^ {p + p) u^Uu +pgnu, (5) 

where p is the energy density and p is the pressure. The Ricci tensor is given in terms 
of the stress tensor by Einstein's equations 



Rixu = Stt [Tfj,^ - - g^yTj , (6) 

and for the perfect fluid we have 

i?^,K'^K" = 47r(p + 3p). (7) 

Thus 9q satisfles 

^ = -4vr(p + 3p) - \el . (8) 
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The conservation law for a perfect fluid can be expressed in terms of the expansion 



as 



26, 27 



p + (p + p) ^ = . 



(9) 



Write the equation of state a.s p = w p and assume w is constant. The density can 
be expressed as 

p{t) = p{0)e-^^+'"^^oHti)dt, _ ^^Q) 
For the case of unperturbed Robertson- Walker spacetime, this is equivalent to 



Pit) = poit) = p(0) 



a(0) 



B. Fluctuations of the expansion 

Now we wish to consider perturbations of the background spacetime produced 
by stress tensor fluctuations. Let 9 = 9o + 9i, p = po + 6p, and p = po + Sp. 
For now, we treat 6p, and 6p as independent variables. Let the Ricci tensor term 
in the Raychaudhuri equation be expressed as a sum of a classical part, given by 

Note that 



Eq. ([7j), and a smaller fluctuating quantum part, denoted by [R^ 
the quantum field responsible for the fluctuations is distinct from the classical perfect 
fluid. Because of the possibility of pressure gradients, the fluid elements will not in 
general move along geodesies and the Raychaudhuri equation acquires an additional 
term on the right hand side (28l . 29| of 



p + p 



(12) 



where is the Laplacian operator in a constant t hypersurface. If we expand the 



Raychaudhuri equation to first order in 9i, 6p, 6p and {R^yU^^u'')q, we find 



d9i 
It 



-47r((5p + 3(5p) 



VHp_ 

Po+Po 



(13) 



This equation may be integrated to find 



6^{t) = -a-'\t) I dt'a\t' 

'to 



An {5 pit') + ?>5p{t')) 



V^5p 
Po+Po 



- {R^^u'^u") 



(14) 



with the initial condition ^i(to) = 0. 

To leading order, we may regard 6p and 6p as perturbations which have some 
source other than quantum stress tensor fluctuations. In most inflationary models, 
these perturbations are driven by quantum fluctuations of the inflaton field, but for 
our purposes they can be treated as being either classical, or at least uncorrected 
with the stress tensor fluctuations. Thus it is convenient to split 9i into a "classical" 
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part 6ic which depends upon Sp and Sp, and a quantum part, 6iq, driven by the stress 
tensor fluctuations. Then 6i = 6ic + 6iq, and 

Oiq{t) = -a-\t) f dt' a\t') {R^uu^u")^ . (15) 

Thus 9iq is given by the fluctuating part of the Ricci tensor, which is in turn given by 
the quantum stress tensor through Eq. ([6]). In this case, we can construct a correlation 
function for the expansion as a integral of the Ricci tensor correlation function: 

K^uaP = {Rfiu{x) Rapix')) - {R^u{x)){Rap{x')) , (16) 

and write 

{e{h)e{t2))-{e{h)){e{t2)) = a-\h)a-\t2) f dta\t) f dt'a\t') K^.^^u^u^u'^u^ 

J to J to 

(17) 

Note that if 6ic is truly classical and non-fluctuating, we do not need to make a 
distinction between 6 and 6iq in the above expression, as the correlation function for 
both is the same: only the fluctuating part contributes to the correlation function. If 
6p and 6p do indeed fluctuate, then we assume that their fluctuations are uncorrelated 
with those of the stress tensor, so their effect would be to add another term in the 
expansion correlation function. This assumption will be discussed in more detail 
later. In writing Eq. (IT7|1 . we are essentially assuming that the ^-fluctuations vanish 
before the t = to hypersurface, which amounts to a sudden switching assumption. 

In this paper, we will restrict our attention to passive metric fluctuations caused 
by conformally invariant quantum fields. In this case, the classical stress tensor 
in Robertson- Walker spacetime is related to that in flat spacetime by a conformal 
transformation: 

T^^ix) = a-\t)Tft^\x). (18) 

The quantum stress tensor operator in curved spacetime has an anomalous trace, 
so curved spacetime expectation values cannot be obtained directly by a conformal 
transformation of the corresponding flat spacetime expectation value. However, the 
contribution of the anomalous trace to the stress tensor operator is a c-number, and 
hence will cancel in a stress tensor correlation function [30]. Thus we can express 
the correlation function for the stress tensor in Robert son- Walker spacetime as the 
conformal transform of the corresponding flat space correlation function, 

C^pix, x) = a-\t) a-\t') Cllt,{x, x') , (19) 

where 

C^.a/3(x,x') = {T^,{x)T^p{x')) - {T^,{x)){T^p{x')) . (20) 

Because the anomalous trace of the stress tensor does not contribute to correlation 
functions, we can use the Einstein equation, Eq. ([S]), to relate the Ricci and stress 
tensor correlation functions: 

K^uap{x, x') = {8nf Cf,yap{x, x') . (21) 
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Here we will assume that the quantum field in Robertson- Walker spacetime is in 
the conformal vacuum state, so the corresponding fiat space correlation function will 
be that for the Minkowski vacuum state. In this case, we have 

K^,^p{x,x')u^u''u''u^ = {MfCut't'{x,x') = {^^fa~\t)a-\t')8, (22) 

where £ is the fiat space vacuum energy density correlation function. 



C. Treatment of singular integrands 

Observable quantities are expressed as integrals of stress tensor correlation func- 
tions, as in Eq. (fT7j) . However, the integrands in these integrals contain singularities 
which appear not to be integrable, specifically higher order poles on the real axis. 
Nonetheless, such integrals can be given an unambiguous, finite value. One approach 
which could be used is dimensional regularization. This approach has been studied 



in fiat spacetime in Ref. [17|, where it was shown that spacetime integrals of stress 
tensor correlation functions such as C^iya/3(x, x') are actually finite in dimensional 
regularization. This means that if we were to evaluate an integral of C^jyQ/3(a;, x') in 
spacetime dimension 4+e, and then take the limit that e — > 0, the result will be finite. 
(This is not true for integrals of correlation functions involving a time-ordered prod- 
uct of stress tensor operators. In this case, the singularity as e is proportional 
to counterterms in the gravitational action quadratic in the curvature.) 

An alternative approach, which is easier to use in practice, is to define the integrals 
by an integration by parts procedure. This is the generalized principal value discussed 



in Ref. [31i, and used in Refs. [12|, [151, Il6|, ll8|. The basic idea is to re-express an 
integral of the form 



(x 



■ dx (23) 



by use of the identity 



1 d^ 
^ = {-lY-\n - 1)! — ln(x - c) , (24) 

and then perform successive integrations by parts to express the original integral as a 
sum of finite boundary terms and an integral whose integrand has only a logarithmic 
singularity. An equivalent approach is to seek an antiderivative, G{x) of the function 
F{x) = f{x)/{x — c)", that is, G'{x) = F{x), and write 

b 

F{x) dx = G{b) - G{a) . (25) 

This result would be trivial if F{x) had no singularities. In the present case, the 
only nontrivial aspect arises from whether the contour of integration goes above or 
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below the higher order pole ai x — c. However, the residue of this pole is pure 
imaginary if f{x) and its derivatives are real. This is the case in the present problem, 
where one needs to integrate a real valued correlation function to get a real answer. 
Consequently, the residue of the pole will not contribute to the final result. 



III. FIRST MODEL: 61-FLUCTUATIONS AND REDSHIFTS AFTER RE- 
HEATING 

In this section, we will examine the fluctuations of the expansion during an in- 
flationary period, and their subsequent effects in creating density fluctuations by 
differential redshifts after the end of inflation. It will be convenient to use the con- 
formal time 7] rather than the proper time t of the comoving observers. The two time 
coordinates are related by 

df] = a~^{t)dt . (26) 

Let the inflationary phase begin at r] = tiq and end at ?7 = 0. During this interval, the 
spacetime will be taken to be deSitter space, which may be represented as a spatially 
flat Robertson- Walker metric with 

= Vo<V<0- (27) 

For ?7 > 0, we take the scale factor to be that of a radiation-dominated universe, for 
which 

a{r])^l + Hr). (28) 

These forms are chosen so that both 0(77) and its first derivative are continuous at 
77 = 0. In terms of comoving time, the scale factors are 

a{t) = e^(*-*«) , t<tR, (29) 

and 

a{t) = ^l + 2H{t- in) , t>tR, (30) 

where t = Ir is the comoving time at which inflation ends. 

The expansion correlation function, both during and after inflation, may be written 

as 

{0{V1)0{V2)) - {0{V1)){0{V2)) = {Snf a-\m)a-\r„) ^ ^ f ^ £{Ar,,r) , 

(31) 

where At] = r) — rj' and r — |x — x'| is the coordinate space separation of the pair of 
points at which 6 is measured. Here we will assume that the ^^-fluctuations vanish 
at the beginning of inflation, r] = tjq. During inflation, there is no classical matter 
present. However, after reheating at 77 = 0, variations in 9 cause the matter in 
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different spatial regions to redshift at different rates, leading to variations in the 
density of the classical matter. Once reheating has occurred, 6p = w 6p, where 
w = 1/3 in our model. Equation ( ITOl) for the energy density can be written as 




(32) 







so that 



- = -(1 + w) f e,{t,) dt, = -(1 +w) 

Po Jtn Jo a{T]i) 



(33) 



The expansion fluctuations lead to density fluctuations given by 




di 



aij]') 



(34) 



In the above expression, the integrals on rji and 772 represent the differential redshift- 
ing, and hence have a lower limit at = 0, the reheating time. The integrals on rj 
and 7]' describe the effects of quantum stress tensor fluctuations on the expansion 9. 
The integration range is the beginning of inflation at 77 = 770 to a time t] = rjs when 
the density variations are measured. We can take rjs to be the time of last scattering, 
when the density fluctuations of the cosmic background radiation were established. 

The expression Eq. (134|1 contains contributions from all length scales, whereas 
observations are sensitive only to a finite range of scales. Thus, in order to compare 
the results of our calculations with observation, we should look at the power spectrum 
Pk{i]s) defined by 



We will first compute the density fluctuations as a function of r, and then Fourier 
transform the result to obtain Pk{rjs). 




(35) 



Let 




Fo + F^ijli) + ^1(^2) + i^2(r?i, V2) , 



(36) 
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where 



Jnn "v/ ) 



Jo (^{V) Jvo '^w ) 

Jvo «(^) ./o aw) 
rfr] rfr/' 



Here Fo describes correlated stress tensor fluctuations entirely within the deSitter 
phase, F2 similarly describes fluctuations entirely in the radiation-dominated phase, 
and Fi describes the correlation of fluctuations between the deSitter and radiation- 
dominated phases. Note that 

F, = {e{o)em-{em{em = {{Aer) (38) 

is the variance of the expansion at the end of inflation. 

To go further, we need the explicit form for the flat space energy density correlation 
function, S{r, Arj). For the case of the electromagnetic fleld, it is 

ii!'(r^ - Arff ■ ^ ' 

For the conformal scalar fleld, 

127r4(r2 - At] 
Let us flrst consider Fq, which may be written as 



r\vo\ rlvol 

Fo= dr]{l + H7]) / d7]' {1 + Hr]')£{Ar],r) . (41) 
Jo Jo 

This integral may be evaluated using algebraic symbol manipulation programs. We 
have used both Mathematica and Maxima with equivalent results. The program 
is asked to flnd antiderivatives of the integrand, which are then evaluated at the 
appropriate limits. (See the discussion in Sec. Ill CI ) The result is rather complicated, 
but simplifles greatly in the limit that i^|?7o| ^ 1- For the electromagnetic fleld case. 



we flnd the asymptotic form 



for if|?7o| ^ 1- For the scalar case in the same limit 

ISvr^r^ 



^o«^ (42) 



(43) 
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The remarkable feature of this result is that it grows with increasing |?7o|, and hence 
depends upon the length of the inflationary period. Note that F2 is independent of 
|?7o|; and Fi is found to go to a finite limit for large It^oI- Thus in this limit, 

F(r?i,r/2) ^Fo. (44) 

Note that 

H\r]o\ ^ e-^(*«-*o) (45) 

is the net expansion factor during inflation, which needs to be greater than about 10^^ 
to solve the horizon problem. Thus the large H\t]q\ approximation is an extremely 
good one. Note that there are no real particles moving along the comoving geodesies 
during inflation. Nonetheless, the expansion 6 at the end of inflation has observable 
consequences. If we consider the reheating to occur very quickly, then we are effec- 
tively matching deSitter spacetime and the radiation dominated Robertson- Walker 
spacetime across the 77 = hypersurface. Any such matching must satisfy the Israel 
junction conditions, that the extrinsic curvature of this hypersurface be continuous. 
(Note that 6 is the trace of the extrinsic curvature tensor of this surface.) This implies 
that 6 must be continuous, even in the case where there are spatial variation in 6. As 
a result, the expansion fluctuations generated by stress tensor fluctuations in the de- 
Sitter phase persist in the radiation-dominated phase and cause density fluctuations. 
Recall that the geodesies in deSitter space with which we are concerned become the 
comoving geodesies in the post-inflationary universe. This choice breaks the deSitter 
invar iance. 

In order to compute the power spectrum of these fluctuations, we must find the 
Fourier transform of which may be done by integration by parts as follows: 

1 1 /"°° 1 

d^x — e"*''''^'' = ^^-j^ I dr — sin(fcr) 



(27r)3 J r6 2n^k Jq r 

^ dr sm{kr) ( —] = — 5- — : / dr ( -] -j-j sin(A;r) 



47r2/c4! J_^ dr^ \r J An'^kAl J _^ \r J dr^ 

r , sin(A;r) k^ 

= / dr = . (46) 

487r2 Jq r 96it ^ ^ 

For the electromagnetic case, this leads to 

, 32H^\rio\^k^ f drii ,2 

+ ,47) 

We may evaluate the integral in the above expression using Eqs. and (1501) to find 

P.iVs) - lAaiVs)] (l + wfil, (48) 

where we have explicitly written the powers of ip, the Planck length. If we are 
interested in the effects of fluctuations within a finite bandwidth, {k, k + Ak), then 
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we can write 

2 



(^^^ j= j d^ke'^-^^PkiVs) ~ /^kk'^e'^-^^Pk. (49) 

For the purpose of a rough estimate, let us take A/c ^ k and e^'^'^^ 1. Then the 
corresponding density perturbation is given by 



Note that we have taken a = 1 at the end of inflation. As a result, l/a{r]s) is the 
redshift factor between reheating and the last scattering surface, 

aiVs) ~ ^ , (51) 
where is the reheating energy scale. We should have 

^1 <10-^ (52) 



p 

which leads to an upper bound on the duration of inflation 



eik^\n{j^) 



If Eji is close to the scale of inflation, then the vacuum energy density during inflation 
is Vo ~ Ej^, and 

^^ = Y4^o-y4^^ (54) 
Because we have chosen a = 1 at the end of inflation, the scale factor today is 



a„ow~10'a(r;.)«103-%, (55) 

lev 

and k is related to the physical wavenumber today, kp, by k = a-^ow kp. Let kp — 
2ti / \ correspond to the typical intergalactic separation today, A ~ 2Mpc, or /cp ~ 
10~^^cm~^. Then we may combine the above relations to write the bound on the 
expansion factor during inflation as 

H|„„|<10-(™). (56) 

ignoring the weak logarithmic dependence upon Er. 
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The dependence found in Eq. (1481) indicates a non-scale invariant spectrum of 
fluctuations which rises at shorter wavelengths. The same dependence upon k was 



found recently in a somewhat different context by Lombardo and Nacir [32|. (See 
Eq. (68) of their paper.) If observational data for kp > 10~^^cm~^ were available, 
then one might be able to obtain tighter bounds on the duration of inflation. The 
smallest scales on which the cosmic microwave background has been observed is about 
5 arcminutes 33|], which corresponds to kp ^ 10~^^cm~^. Similarly, the role of higher 



values of fcp in large scale structure formation is unclear because of nonlinear classical 
effects. 

The effect we are considering depends upon transplanckian modes in the sense that 
the modes of the quantized electromagnetic or scalar fields which give the dominant 
contribution have wavelengths much shorter than the Planck length. Let Lj be a given 
proper length at the beginning of inflation, and let Lf he the corresponding scale 
today stretched by the cosmological expansion. These two scales are approximately 
related by 

^/ = 10' (^) ^l^ol^- (57) 

as H ItjqI is the expansion during inflation, Eji/leV is that between reheating and 
last scattering, and there has been an additional expansion by a factor of about 10^ 
after last scattering. If we take Lf ^ lO^^'cm, Eji ^ lO^^GeV, and H |?7o| ~ lO''^, then 
Li ^ 10~^^£p. It is well known that transplanckian modes play a crucial role in the 
conventional approach to black hole evaporation [3^]. It is possible to obtain black 
hole evaporation without transplanckian modes, but only at the price of introducing 



a Lorentz non-invariant dispersion relation [35 . 



IV. SECOND MODEL: SINGLE FIELD SLOW-ROLL INFLATION 

Let be an inflaton field which obeys the equation of motion in a Robertson- 
Walker metric 

(58) 

where V{(j)) is a relatively flat potential, and is the flat space Laplacian operator. 
During a slow-roll phase, the second time derivative of (p is assumed to be small. If, 
in addition, the spatial gradient terms are small, then 

3-^^V'{(j)). (59) 
a 

Note that = 3d/a = 3H is the expansion of unperturbed deSitter spacetime. 
Now we wish to generalize this description to include small spatial variations of the 
expansion. If the vorticity of the comoving geodesies vanishes and the shear remains 
small, then the spacetime metric can be written as 

ds'^ = -dt^ + a^{t, x){dx'^ + dy^ + dz^) . (60) 
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The expansion of the comoving geodesies, those with four-velocity u'^ = in these 
coordinates, is again 

^ = < = 3-, (61) 

and the shear and vorticity vanish. The equation of motion for an inflaton field with 
self-coupling V{(f)) now becomes 

□^=J_V.(„V^)-l|(„»f) = rW. (62) 

Let each of the quantities a, 0, and 9 consist of a homogeneous part and a small 
inhomogeneous perturbation: 

a = ao(t) + ai()f:, x), = 0o(^) + x), and 6 = Ooit) + 9i{t,:s.) . (63) 

We make the slow- roll approximation for the homogeneous part, (poit), which satisfies 

^o^-q^. (64) 

where = 3do/ao- We retain the second time derivative of (j)i{t,x), which satisfies 

01 + ^0 01 + ^1 00 - ^ V201 = -V"{<Po) 01 , (65) 

where we have expanded Eq. (!62l) to first order in all of the inhomogeneous pertur- 
bations and used 

y (00 + 0i) ^ V'(0o) + V"(0o) 01 . (66) 

Let us consider the case where the potential is approximately linear during the period 
of interest, so we may set V"{(j)o) ~ 0. If we use Eq. f lM|) and 6q = 3H, we may write 
Eq. §E> as 

1 -V2. V'{<Po) 



01 + 3/701 - - V>i = — ^ 01 . (67) 
Let us Fourier transform this equation and define 

SMt) = j ci3A;e^''-^^0i(t,x), (68) 



and 



SOkit) = j d^ke'^-^''eiit,x) . (69) 
We also change from comoving time t to conformal time rj. Then 6(f)k{v) satisfies 

-W + + k5cp, = ^a, 50, . (70) 
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Let G{ri, rj') be a retarded Green's function for this equation which satisfies 

"""^ +2Hao^ + k'G = 5{rj-r^'), (71) 



dr]"^ dr] 
and G{ri, 77') = if 77 < r]'. It can be expressed as 

Here ipi and are two hnearly independent solutions of Eq. (17X1) without the delta- 
function source term, and W is their Wronskian. Here t]^ and ?7< are, respectively, 
the greater and lesser of rj and t]'. These solutions may be taken to be 

Mv) = Iv^h\v\Ih^:^\v) (73) 

z 2 

and 

Mv) = lV^H\r,\lHf\7i), (74) 
z 2 

where iyl^"* and H^^^ are Hankel functions. The Wronskian becomes 



3 aiiu. -I J 3 

2 2 



If we set 



H'bi('/)*'2('/)l = ^. (75) 



F = ^<.?«., (76) 



then we can write the solution of Eq. (!70l) as 

6MV) = f dv' F{r^') G{r^, i) = P d^ F{i) G{r^, V) , (77) 



where we use the fact that G{ri,ri') = for rj' > rj. In the expression for G{ri,ri'), we 
may now set f]^ = rj and ?7< = 77'. Let us now split the integration range into two 
parts as 

Hk{ri) = r dri' F{r]') G{r], r]') + /" dr]' F{r]') G{r], r]') . (78) 



vo 



Here r/^ = — 1/fc is the horizon-crossing time for mode k in conformal time. Note 
that in the first integral in the above expression, we have klrjl < 1 < k\ri'\, and in the 
second integral we have k\r]\ < k\r]'\ < 1. If we assume that k\r]\ <^ 1 -C k\r]'\, and 
use the limiting forms for the Hankel functions, we find 



G=»^. (79) 
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Similarly, if we assume k\r]\ <^ k\ri'\ <^ 1, then we have 



G 



^0) 



Let us now make the approximation that we may use Eq. (!79|) in the first integral in 
Eq. (ITHIl . and Eq. (IHOll in the second. The result is 



9H^ 



'13 



n 1 ' 



56, 



^1) 



where we have used k = l/|?7c| and our assumption that Vq = V'{(f)o) is approximately 
constant. 

We now consider the effects of expansion fluctuations upon the evolution of the 
inflaton field. The variance of S(j)k{Ti) is 



Vc 



3{Hr]c) / rfr/i ao(?7i) cos(r7i/?7c) - / drji 



vo 



Vc 



Vc 



dr]2al{r]2)cos{T]2/T]c) - / rfr/; 



Vo 



Vc 



X ao(r/i)ao(r72) [{69^ 50^) - {S9k{vi)){5ek{V2))] ■ 



(82) 



The correlation function in the above integral is just the Fourier transform of the 
coordinate space expansion correlation function given in Eq. (13T|) . It is convenient to 
evaluate the integrals in coordinate space and to write variance of S(j){ri) as 



((A0)^) 



X 



\9H\ 



Vc 



dr]ial{r]i)cos{r]i/r]c) - / dr] 
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where 
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vo 



V2 



dvao\v) / dr]' aQ\r]')£{A7],r) 



(83) 



vo 



Here we have evaluated the variance of at the end of inflation, 77 = 0. This 
quantity is directly related to the density perturbations in the post-inflationary era. 



This issue has been discussed extensively by previous authors [36|, l37|, l38|, l39|, |40 



m a 



context where A0 arises from the intrinsic quantum fluctuations of the inflaton field. 
However, the relation between A0 and the density perturbation is the same in the 
present context. The basic idea is that spatial variations in cause different regions 
to reheat at different times, with a typical time variation of order At ~ A0/0o- After 
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reheating, the energy density decreases as p = pr (t/j/t)^, where pr is the density at 
reheating. This leads to density variations whose magnitude is of order 

6p 2At 2A(j) 6HA(f) 

— ^ ^ ^ — - . (85) 

P tR tR(f)0 tRVQ 

In the last step, we have used Eq. (!64l) with = V(0o)- Thus we have 

2 




6H 



((A0)^), (86) 



with ((A0)^) given by Eq. ( !83|) . The physical picture for the conversion of 0- 
fluctuations into density variations described above was first given by Guth and 



Pi |37[. It is possible to give a more rigorous, gauge- invariant discussion [29|, 141 
However, for a single inflaton model, the results are essentially equivalent to Eq. flHB]) . 

Consider first the case of an electromagnetic field. It is convenient to express 
Eq. ([39D as 



S - 1 



4807r4 



(r2 + 3Ar/' 
ds^ {s — At]' 



2\2 



2^ 



^7) 



and to interchange the order of the s-differentiations and the integrations in Eq. 
In the limit of large \riQ\, the result is 



I ^ JL^ , (88) 

which is independent of rji and ri2 to leading order. Now the integrations on these 
variables in Eq. (183!) can be performed in terms of the cosine integral function, ci{x), 
with the result 

where ci{l) ~ 0.337. Here we assume that H\t]c\ ^ 1. The corresponding expression 
for the conformal scalar field is smaller by a factor of 1/3. In this model, the 
fluctuations grow at the same rate with increasing |?7o| as in the case of the model in 
Sect, mil This arises despite the fact that A0fc for an individual mode only begins 
to grow after the mode leaves the horizon. The growth with increasing |?7o| comes 
from the fact that the expansion fluctuations grow from the beginning of inflation, 
as described by Eq. flSTl) . If it were not for this growth, the factor of ip would 
tend to make the effects of stress tensor fluctuations small compared to those of the 
intrinsic quantum fluctuations of (j). The corresponding power spectrum of density 
perturbations is 

M tt2 2 4 U3 pi tt2 2 
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where we set |?7c| ~ 1/fc only at the last step. If we make the same estimates as were 
made in Eq. ( !50l) . we find 



\ P J rms 



tR 



Note that both here and in the model of Sect. IIIIl we find a spectrum which is not 
scale invariant. Here {Sp/p)rms k. If we again set kp ~ lO^^^cm^^ and assume 

that Eji ^ Vq^ , we obtain the following upper bound on the duration of infiation: 

^|,„|<10«fl2^)'. (92, 



E 



R 



This is considerably more restrictive than Eq. (156!) . but is still compatible with ade- 
quate infiation to solve the horizon problem. 



V. SUMMARY AND DISCUSSION 

We have analyzed the effects of stress tensor fluctuations of conformally invari- 
ant quantum flelds in deSitter spacetime. One unexpected result of this analysis 
is that expansion fluctuations grow during a deSitter phase. This might be inter- 
preted as due to the background spacetime altering the anti-correlated fluctuations. 
In Minkowski spacetime, quantum fluctuations tend to have strong correlations and 
anti-correlations. If one were to evaluate the variance of the expansion in Minkowski 
spacetime, that is, compute Eq. (!3T!) with a = 1, the result would be independent of 
?7o in the limit that Ir^ol oo. This is closely related to the ant icor relations found 
in sampled energy density measurements in flat spacetime. Similar anticorrelations 
are present when a charged particle or a mirror is coupled to vacuum fluctuations in 
flat spacetime, causing the mean squared velocity to approach a constant even in the 



absence of dissipation [42|, |43|, |4^. It is the presence of nonconstant 0(77) functions 
in the integrand of Eq. fl3ip which upsets the flat spacetime cancellations and leads 
to a result which grows with increasing Ir^ol- 

It is this growth which allowed us to infer the constraints, Eqs. fl56l) and fl92l) on 
the duration of the inflationary phase. These constraints may come as a surprize, as 
one usually expects inflation to erase the memory of the past history of the universe. 
It is certainly true that the exponential expansion quickly suppresses classical pertur- 
bations. However, the effect discussed here amounts to a type of quantum instability 
of deSitter spacetime: the cumulative effects of passive metric fluctuations eventu- 
ally lead to a spatially inhomogeneous spacetime. The direct effect on the spacetime 
geometry grows rather slowly, as reflected in the constraint Eq. (jSSD- However, when 
the ^-fluctuations couple to the inflaton fleld, the result can be a stronger constraint, 
Eq. (p2!) . Both of these constraints are consistent with adequate inflation to solve 
the horizon problem. The possibility of effects which grow during inflation and react 
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against the expansion has been discussed by several authors as a possible solution 
for the cosmological constant problem. Amon g th e effects considered are the growth 
of long waveleiigth classical perturbations (isl . 4^ and backreaction due to quantum 
gravity effects |47(]. The effect discussed in the present paper is distinct from either of 
these effects, especially as it produces an increasingly inhomogeneous universe rather 
than a backreaction against the cosmological constant. 

Here is it worthwhile revisiting three of the assumptions which we made in our 
analysis. One was the assumption that the 5p and 5p terms in Eq. f|T^ are uncorre- 
cted with the term produced by stress tensor fluctuations. This assumption should 
hold so long as the dominant source of density perturbations is other than the stress 
tensor fluctuations. Of course, if inflation were to last sufficiently long, this would no 
longer be true; 5p and 5p would be predominately due to these fluctuations. How- 
ever, because the contributions of stress tensor fluctuations are highly non-Gaussian 
and non-scale invariant, they can give at best a very small contribution to the total 
primordial fluctuation spectrum in our universe. Hence our assumption seems to be 
justifled. Another assumption which we have made is ignoring the shear term in the 
Raychaudhuri equation. The comoving geodesies in deSitter spacetime are certainly 
shear-free, so it is reasonable to start in a state where cr^jy = 0. Shear could de- 
velop only if there were large fluctuations of the Weyl tensor. This seems unlikely, 
but Weyl tensor fluctuations from passive metric fluctuations need to be better un- 
derstood. The third assumption is the sudden switching assumption flrst made in 
writing Eq. fll7l) . This amounts to assuming that the expansion fluctuations vanish 
before the t = to hypersurface, and that the effects of the stress tensor fluctuations 
appear suddenly after that time. In future work, we plan to examine more general 
initial conditions to test the dependence of our results upon the initial conditions. 

A particularly interesting possibility is that inflation lasted for a time only slightly 
less than the constraints derived above. In this case, one would predict a small, but 
potentially observable effect from the quantum stress tensor fluctuations. This effect 
is expected to manifest itself in a non-Gaussian and non-scale invariant component 
in the density perturbations. It would also offer a possible probe of transplanckian 
physics. 



Acknowledgments 

We have benefltted from discussions with numerous colleagues. In particular, 
L.H.F. would like to thank the participants of the 11th Peyresq workshop for valu- 
able comments. This work was supported in part by the National Science Council 
of Taiwan under Grant NSC95-2112-M-001-052-MY3, and by the National Science 



18 



Foundation under Grant PHY-0555754. 



[1] 
[2] 
[3] 
[4] 
[5] 
[6] 
[7] 

[8] 
[9] 

[10] 

[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 



L.H. Ford, Ann. Phys (NY) 144, 238 (1982). 

Chung-I Kuo and L.H. Ford, Phys. Rev. D 47, 4510 (1993), |gr-qc/9304008'. 



B.L. Hu and K. Shiokawa, Phys. Rev. D 57, 3474 (1998), |g r-qc/9708023| 
B.L. Hu, Int .J. Theor. Phys. 38, 2987 (1999), li?^q^/9902064j 



N.G. Phillips and B.L. Hu, Phys. Rev. D 63, 104001 (2001), ^r^qc/OOlOOlg 



B.L. Hu and E. Verdaguer, Class. Quant. Grav. 20, Rl (2003), gr-qc/0211090 



E. Calzetta and B.L. Hu, Phys. Rev. D 49, 6636 (1994), |gF^c/9312036t Phys. Rev. D 
52, 6770 (1995), gr-qc//9505046. 

E. Calzetta, A. Campos, and E. Verdaguer, Phys. Rev. D 56, 2163 (1997), 



gr-qc/9704010 



R. Martin and E. Verdaguer, Phys. Lett. B 465, 113 (1999), gr-qc/9811070 Phys. Rev. 



D 60, 084008 (1999), |gr-qc/9904021|; Phys. Rev. D 61, 124024 (2000), gr-qc/0001098 



B.L. Hu and E. Verdaguer, Living Rev. Rel. 7, 3 (2004), .gr-qc/0307032 



J.W. Moffat, Phys. Rev. D 56, 6264 (1997), gr-qc/9610067 , 



J. Borgman and L.H. Ford, Phys. Rev. D 70 064032 (2004), lgr-qc/0307043 



J. Borgman and L.H. Ford, Phys. Rev. D 70, 127701 (2004), !gr-qc/0406066 



C.-H. Wu and L.H. Ford, Phys. Rev. D 60, 104013 (1999), lgr-qc/9905012 



C.-H. Wu and L.H. Ford, Phys. Rev. D 64, 045010 (2001), quant-ph/0012144 



L.H. Ford and C.-H. Wu, Int. J. Theor. Phys. 42, 15 (2003), gr-qc/0102063 , 



L.H. Ford and R.P. Woodard, Class. Quant. Grav. 22, 1637 (2005), gr-qc/0411003 



L.H. Ford and T.A. Roman, Phys. Rev. D 72, 105010 (2005), gr-qc/0506026 



C.-H. Wu, K.-W. Ng, W. Lee, D.-S. Lee, and Y.-Y. Charng, JCAP 2, 6 (2007), 



|astro-ph/0604292 , 

A.H. Guth, Phys. Rev. D 23, 347 (1981). 

See, for example, A.R. Liddle and D.H. Lyth, Cosmological inflation and large-scale 



19 



structure (Cambridge University Press, 2000). 
[22] D. N. Spergel et al, | astro-ph/06 03449} 

[23] A. Borde, A.H. Guth, and A. Vilenkin, Phys. Rev. Lett. 90, 151301 (2003), 



gr-qc/0110012 



[24] L.H. Ford, quant-ph/0601112 



[25] C.J. Fewster, L.H. Ford, and T.R. Roman, manuscript in preparation. 

[26] S.W. Hawking, Astrophys. J. 145, 544 (1966). 

[27] D.W. Olson, Phys. Rev. D 14, 327 (1976). 

[28] D.H. Lyth and E.D. Steward, Astrophys. J. 361, 343 (1990). 



[29] A.D. Liddle and D.H. Lyth, Phys. Rep. 231, 1 (1993), |astro-ph/9303019 



[30] N.D. Birrell and P.C.W. Davies, Quantum fields in curved space (Cambridge University 
Press, 1982). 

[31] K.T.R. Davies and R.W. Davies, Can. J. Phys. 67, 759 (1989); K.T.R. Davies, R.W. 
Davies, and CD. White, J. Math. Phys. 31, 1356 (1990). 



[32] F.C. Lombardo and D. Lopez Nacir, Phys. Rev. D 72, 063506 (2005), [gr^/0506051 
[33] C.L. Kuo et al, Astrophys. J. 600, 32 (2004), | astro-ph/0212289[ 
[34] S.W. Hawking, Commun. Math. Phys. 43, 199 (1975). 



[35] S. Corley and T. Jacobson, Phys. Rev. D 54, 1568 (1996), [hep^h/9601073 

[36] V. Mukhanov and G. Chibisov, JETP Lett. 33, 532 (1981). 

[37] A.H. Guth and S.-Y. Pi, Phys. Rev. Lett. 49, 1110 (1982). 

[38] S.W. Hawking, Phys. Lett. B 115, 295 (1982). 

[39] A.A. Starobinsky, Phys. Lett. B 117, 175 (1982). 

[40] J.M. Bardeen, P.J. Steinhardt, and M.S. Turner, Phys. Rev. D 28, 679 (1983). 

[41] V.F. Mukhanov, R.H. Brandenberger, and H. Feldman, Phys. Rep. 215, 203 (1992). 



[42] C.-H. Wu and D.-S. Lee, Phys. Lett. A 318, 303 (2003), |gr-qc/0304034, 



[43] C.-H. Wu and D.-S. Lee, Phys. Rev. D 71 125005 (2005), quant-ph/0501127 



[44] H. Yu and L.H. Ford, Phys. Rev. D 70 065009 (2004), [qn^^ph70406122] 



[45] G. Geshnizjani and R. Brandenberger, Phys. Rev. D 66, 123507 (2002), gr-qc/0204074 



20 



[46] R. Brandenberger and A. Mazumdar, JCAP 0408, 015 (2004), |hep-th/0402205 



[47] N. C. Tsamis and R. R Woodard, Annals Rhys. 253, 1 (1997), hep-ph/9602316, 



21 



